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Abstract 

We present results for ir, p and proton and A masses from our recently built Pentium 
cluster. The previous results for quenched Wilson fermions by MILC and GF11 collaborations 
are compared at (3 = 6/g 2 —5.7 and 5.85 with the same parameters on the 8 3 x 32 lattice. 
New data are shown at (3—6.25 on the 12 3 x 36 and 16 3 x 32 lattices. Such a larger (3 value is 
useful for extrapolating the lattice results to the continuum limit. The smearing technique is 
systematically investigated and shown to greatly improve the spectrum data. 
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1 Introduction 



The computational physics group at the Zhongshan University is in a period of rapid develop- 
ment. The group's interests]^ [| cover such topics as lattice study of hadron spectroscopy []|, M, 
glueball decay and mixing, QCD at finite density ||, quantum instontons[0] and quantum chaos|§] 
using the quantum action||. Most of these topics can be studied through Monte Carlo simula- 
tion, but can be quite costly in terms of computing power. In order to do large scale numerical 
investigations of these topics, we built a high performance parallel computer]!!], 11] using PC 
components. 

QCD has been accepted to be the most successful gauge theory of strongly interacting particles. 
Calculation of hadron spectroscopy remains to be an important task of non-perturbative studies 
of QCD using lattice methods. This paper is one of our first steps || |l0| in the direction of 
studying hadron properties. Our reasons for performing simulations with Wilson valence quarks 
are twofold: First, we are interested in the analysis of hadronic spectroscopy for quenched Wilson 
fermions; Second, we are interested in exploring the performance of our cluster in actual lattice 
QCD simulations. We decided to do a spectrum calculation on the 8 3 x 32 lattice at (3 = 5.7 and 
5.85, and on the 12 3 x 36 and 16 3 x 32 lattices at f3 = 6.25. The results at smaller (3 values are 
used to compare with those in the literature, while those at a larger /3 value, corresponding to 
smaller lattice spacing a, are useful for extrapolating the lattice results to the continuum limit. 
The smearing technique is employed to improve the spectrum data. We hope our data will be an 
important addition to the lattice study of QCD spectrum. 

The remainder of this paper is organized as follows. In Sect. ^, the basic ideas of lattice QCD 
with Wilson quarks are given. In Sect. £3, we describe the simulation parameters as well as some 
basic features of our cluster. In Sect. |3], we present new calculation of some light hadron masses. 
Finally, conclusions and outlook are given in Sect. §. 

2 Basic ideas of lattice QCD 

Our starting point is the Wilson action |[^] 

s=-~E Tr (^ + u t - 2 ) + E i>{x)M xy m, (i) 

p x,y 

where f3 = 6/g 2 , and U p is the ordered product of gauge link variables U^(x) = e 9 ' x ' 

around an elementary plaquette. a is the lattice spacing, fx is the unit vector and M is the 
fermionic matrix: 

4 

M xy = 5 xy - k E [(! ~ 7^)U l _ l (x)5 Xt y^f l + (l + ^ f^U^x - fi)5 x ,y+ji ■ (2) 



The fermion field ip on the lattice is related to the continuum one tp cont by ip = r ip cont ^a^'/ (2/t) 
with k = l/(2ma + 8). Then all the physical quantities are calculable through Monte Carlo 
(MC) simulations with importance sampling. Fermion fields must be integrated out before the 
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simulations, leading to the vacuum expectation value for an operator F 

J[dU]F([U])e- s ^(detM) N f p 1 ^ 
W J[dC/]e-^(^)(detM)^' ^iVconfig^ 1 J " lJ 

Here Nf is the number of flavors, and F is the corresponding operator after Wick contraction 
of the fermion fields and the summation is over the gluonic configurations, C, drawn from the 
Boltzmann distribution. In quenched approximation, det(M) = 1. 
The correlation functions of a hadron is: 

C h (t) = J2(Ol(x,t)O h (6,0)), (4) 

X 

where Oh(x,t) is a hadron operator. For sufficiently large values of t and the lattice time period 
T, the correlation function is expected to approach the asymptotic form: 

C h (t) -> Z h [exp(-m h at) + exp(m h at - m h aT)]. (5) 

Fitting the large t behavior of the correlation function according to Eq. (||), the hadron mass mh 
is obtained. 



3 Simulations 

The quenched gauge field configurations and quark propagators were obtained using the recently 



built ZSU cluster] 1C, 11|. The full machine has 10 nodes, and each node consists of dual Pentium 
CPUs. The total internal memory is 1.28GB and the sustained speed is around 2 Gflops. An 
upgrade is planned. The cluster runs on the Linux operating system. To operate the cluster as 
a parallel computer, the programmer must design the algorithm so that it appropriately divides 
the task among the individual processors. We use MPI (Message Passing Interface), one of the 
most popular message passing standards. 

We updated the pure SU(3) gauge fields with Cabibbo-Marinari quasi-heat bath algorithm [p^], 
each iteration followed by 4 over-relaxation sweeps. The simulation parameters are listed in Tab. 
|. The dist ance between two nearest stored configurations is 100. The auto-correlation time was 
computed to make sure that these configurations are independent. 

The u quark and d quark are assumed to be degenerate. The quark propagators are calcu- 
lated in the Coulomb gauge using the independent configurations mentioned above, and conju- 
gate gradient for inversion of the Dirac matrix with preconditioning via ILU decomposition by 



checkerboards |l4fl . 



At j3 = 6.25, we computed the tt, p and proton and A masses on the 12 3 x 36 and 16 3 x 32 
lattices at four values of hopping parameter: k = 0.1480, 0.1486, 0.1492, and 0.1498 with point 



and smeared[17, 18] sources. In order to compare our results with those by MILC and GF11, the 
quenched simulations were also performed at (3 = 5.7 and (3 = 5.85 on the 8 3 x 32 lattice. We 
repeated the quenched simulations of the MILC and GF11 collaborations, using the same set of 
k values |l6|, 17, O] but on a 8 3 x 32 lattice. Detailed data will be presented elsewhere [10]. 
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4 Light hadron masses 



To extract masses from the hadron propagators, we must average the correlation function of the 
hadron over the ensemble of gauge configurations, and use a fitting routine to evaluate the hadron 
masses mn- We determined hadron masses by fitting our data under the assumption that there 
is a single particle in each channel [p~C| ] . 

Point source means a delta function, and smeared source means a spread-out distribution (an 
approximation to the actual wave- function of the quantum state). For example, the simplest 
operator for a meson is just Oh{x) = q(x)q(x), i.e. the product of quark and anti-quark fields 
at a single point. A disadvantage of this point source, is that this operator creates not only the 
lightest meson, but all possible excited states too. To write down an operator which creates more 
of the single state, one must "smear" the operator out, e.g. J2y Q{%)f(x ~ Vjliv) where f(x) is 
some smooth function. Here we choose 

/(x) = iVexp(-|x| 2 /r 2 ), (6) 

with N a normalization factor. The size of the smeared operator should generally be comparable 
to the size of the hadron created. There is no automatic procedure JL7| ] for tuning the smearing 
parameter tq. One simply has to experiment with a couple of choices. 

The effective meson mass for ir and p as a function of time t at (3 = 6.25 and on the 16 3 x 32 
lattice is depicted in Figs, [l] and [2] respectively. As one sees, the plateau from which one can 
estimate the mass, is very narrow for point source, due to the reason mentioned above. When the 
smearing source is used, the width of the plateau changes with the smearing parameter r$. We 
tried many values of ro and found that when ro > 16, the effective mass is almost independent 
on ro and the optimal value for ro is 26, where we observe the widest plateau. 



The 7r, p, proton and A masses [10] at (3 = 5.7 and on 8 3 x 32 are consistent with the previous 
GF11 data g§ and those for the 

mesons at (5 = 5.85 and on 8 3 x 32 are comparable with those 
in 18] on the 16 3 x 32 lattice. This means the finite size effects are very small. At (3 = 5.85, 
we also see that our results for the proton and A masses |To| are bigger than the MILC data and 
consistent with the finite size behavior analysis in [|l7j . More detailed results will be given in 0] . 
In Fig. ||, we show the Edinburgh plot, {m^ /m p ) vs. (m^/rrip) mass ratios. 

In Figs. [I], and ||, we show respectively (m^a) 2 , m p a, m p a and m&a as a function of 1/k at 
P = 6.25 and on the 16 3 x 32 lattice. Assuming that (m n a) 2 is linear in 1/k, we can compute the 
critical coupling k c at which the pion becomes massless. We extrapolate the data using: 

(m^a) 2 =a(---). (7) 

\K K C J 

The results for k c , and m p a, mpa and m^a at k c on the 12 3 x 36 and 16 3 x 32 are given in Tab.|2[ 
where the experimental value for the p mass has been used as an input. The finite size effects are 



consistent with [17]. Also, the mass ratios mp/m p and m/±/m p seem to be closer the experimental 
values, 1.222 and 1.604 on the larger lattice. To compare mp and ttia with experiment, we need 
to do simulation on a larger lattice volume and carefully study the lattice spacing errors. In this 



aspect, it might be more efficient to use the improved action and some progress |20], [21|] has been 
reported by some lattice groups in China. 
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5 Conclusions and Outlook 



In this paper, we have presented new data on hadron masses in QCD with Wilson valence quarks 
at /3 = 6.25 on the 12 3 x 36 and 16 3 x 32 lattices. Our results at (3 = 5.7 are consistent and those 
at (3 = 5.85 are comparable with the results in the literature. We have also made a more careful 
and systematic study of the smearing method. Such large scale simulations had usually required 
supercomputing resources [jlj, 18], but now they were all done on our recently built cluster for 



parallel computer. In the task of lattice QCD simulations, we are confident that ZSU's Pentium 
cluster can provide a very flexible and extremely economical computing solution, which fits the 
demands and budget of a developing lattice field theory group. We are going to use the machine 
to produce more useful results of non-perturbative physics. 
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volume 


P 


warmup 


stored configs. 


8 3 x 32 


5.7 


200 


200 


8 3 x 32 


5.85 


200 


200 


12 3 x 36 


6.25 


200 


200 


16 a x 32 


6.25 


600 


600 



Table 1: Simulation parameters 



Lattice Volume 


12 3 x 36 


16 3 x 32 


K c 


0.1531(7) 


0.1539(11) 


a(K c ) 


0.505(9) GeV" 1 


0.520(5) GeV- 1 


m 7r a(K c ) 


0.0 


0.0 


m p a(K c ) 


0.388(6) 


0.399(8) 


mpa(K c ) 


0.621(8) 


0.611(5) 




0.770(3) 


0.763(8) 


mp/m p 


1.598(4) 


1.529(8) 




1.981(9) 


1.908(9) 



Table 2: Results of fitted k, lattice spacing, and light hadron masses in the chiral limit at /3=6.25. 
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Figure 1: Pion effective mass fits to the correlation function at (3 = 6.25 and k = 0.1480 and 
on the 16 3 x 32 lattice. Data for the point source, smearing source for ro = 8, and ro = 26 are 
labelled by circles, diamonds, and squares respectively (from top to bottom). 
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Figure 2: The meaning of the symbols are the same as Fig. |l|, but the rho meson. 
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Figure 3: Edinburgh plot. Data are: simulations at (3 = 5.7 (star), j3 = 5.85 (filled circle) on the 
8 3 x 32 lattice, = 6.25 on the 12 3 x 36 lattice (diamond) and [3 = 6.25 on the 16 3 x 32 lattice 
(filled square). The open square shows the real- world point. 
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Figure 4: Pion mass squared as a function of 1/k for f3 = 6.25 on the 16 3 x 32 lattice. In this 
and the following figures, the error bars reflect the sum of statistical and systematic errors. 
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Figure 5: Rho, proton and Delta masses as a function of 1/k for (3 = 6.25 on the 16 3 x 32 lattice. 
Results in the chiral limit (on the vertical dot-dashed line) are also shown. 
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